Structural Aspects of Magnetic Coupling in CaV409 
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The strong corrugation of the V40g layer in the spin gap system CaV409 is examined for its 
impact on the exchange coupling constants between the spin i V ions. Local spin density (LSD) 
calculations show that the V spin occupies a d^2_y2 orbital {x and y are the V-V directions). The 
Kanamori-Goodenough superexchange rules, and the degeneracy of ferromagnetic and antiferromag- 
netic alignments within LSD, suggest that 2nd neighbor couplings dominate over nearest neighbor, 
resulting in two coupled S=i metaplaquette systems. 



From its vanishing susceptibility x{T ^ 0) [|l| and re- 
lated NMR properties |Q CaV409 is known to enter a 
quantum disordered phase with a spin gap. This be- 
havior has stimulated theoretical investigation of the ex- 
change couplings between the S=^ spins on the V lattice, 
[^-^ using a Heisenberg model with nearest (J„„) and 
■ next-nearest ( J„„n) exchange couplings. The ^-depleted 
, lattice (described below) has been viewed as an array 
' of square "plaquettes" of V ions tending toward singlet 
C . formation, g]. Isolated plaquettes have a singlet ground 
state; structural chemistry however suggests intra- and 
inter-plaquette V-V coupling (J and J') should be very 
similar, so the limit of isolated plaquettes is not realistic. 

Depletion alone does not destroy Neel order, |^ and 
although competing interactions lead to x(T)^0, they 
do not account quantitatively for the x(T) data. [^] 
Yet the "plaquette phase" of incipient singlets provides 
an attractive framework for accounting for the lack of 
magnetic ordering and the spin gap behavior. In this pa- 
per we take into account all complications of the crystal 
structure, and find that the orbital character of the oc- 
cupied V d spin orbital is different than previously antic- 
ipated. For this spin orbital, nnn exchange cou- 
pling exceeds nn coupling. The resulting coupling leads 
to a picture of two coupled metaplaquette systems that 
may preserve the tendency toward local singlet formation 
while enhancing frustration of ordering. 

The crystal structure of CaV409 jl^] has been ideal- 
ized in most previous theoretical treatments to the con- 
sideration of VO layers with periodic V vacancies. The 
actual structure is much more interesting. The space 
group is simple tetragonal P4/n (#85 in the International 
Tables), with a = 8.333A, c — S.OOSA, and two formula 
units per cell. The low crystal symmetry is reflected in 
the fact that ten of the 15 internal structural parameters 
(for five sites) are not fixed by symmetry. Conceptually, 
one may start from a VO square lattice with cations and 
anions arranged as on a checkerboard. Each V ion along 
a V-0 axis has an apical oxygen above (say) the layer, 
with V ions along neighboring V-0 lines having their 
apical oxygen below the plane. V ions are removed in a 
pattern corresponding to a left-moving (or right-moving, 
giving the chiral partner) knight on a chessboard that lies 
at a 45° angle derived only from the V ions: from one V 



vacancy, go two V ions along a line of V ions, and one to 
the left to locate a neighboring vacancy. The resulting 
pattern, shown in Fig. 1, is a ^/b x ^/5 enlargement with 
respect to the V sublattice, which is itself x \/2 larger 
than the VO square lattice. Ca ions arrange themselves 
either above or below the V vacancies. 




FIG. 1. Top view of the V4O9 layer, with Ca and 02 atoms 
not displayed. A plaquette of V ions is labelled by 'V. Two 
metaplaquettes, formed by the V ions above the layer, are 
shown by the solid lines. 

Superimposed on this superstructure is a severe cor- 
rugation of the plane, with V ions displaced alternately 
±0.625 A perpendicular to the plane in the direction of 
their apical O ion. The entire structure is shown in Fig. 
2. V ions, as well as apical ions (02), remain equivalent 
by symmetry, but there are two other O sites. 01 sites 
are near (±0.12 A) the plane and are coordinated with 
three V ions, two on one plaquette and the third on a 
neighboring plaquette, all at 1.95-1.96 A. 03 sites remain 
in the plane at positions at the plaquette center, and are 
four-fold coordinated with V ions at 2.04 A. The com- 
pound can be written descriptively as [CaV40l402403]2. 
The V, 01, and 02 sites have no symmetry. The V ion 
lies within a tilted square pyramid of O ions ('square' of 
OI3O3, with 02 apex), not far from the center of mass 
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of the five O ions (Fig. 3). Each V4 plaquette is boat- 
shaped, with two corner V ions up and the other two 
down. Due to the corrugation, the V-V-V angle that 
would be 180° for a plane is reduced to 130°. 




FIG. 2. Structure of CaV409, from Bouloux and Galy, il- 
lustrating the severe corrugation of the V-0 layer. Unit cells 
in two successive layers are shown. 
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the O 2p states to be centered 6-7 eV below the center of 
the V 3c? bands, and separated from them by a gap of 2 
eV. The lowest Ca-derived band bands arise from the 3d 
states, but lie above the V d bands and are inactive. The 
Fermi level Fip lies low in V bands, corresponding to a 

configuration. Thus the conventional ionic picture is 
well respected in this compound. 

To interrogate magnetic interactions, ferromagnetic 
(FM) and antiferromagnetic (AFM) alignments of the V 
moments were studied. The energy gain from polariza- 
tion is 0.12 eV/V ion. For FM ordering, an insulating 
state is obtained, as shown in Fig. 4, so the moment ob- 
tained is precisely 1 per V. An insulating result was 
not anticipated from the non-polarized calculation: there 
was no gap in the d bands, so a rigid Stoner splitting of 
majority and minority bands would leave a metallic re- 
sult. The lowest lying V d orbital in the majority bands, 
which we show below to be the d^2_y2 orbital, upon po- 
larization becomes separated from the remaining 4 d or- 
bitals, due to an exchange splitting A^x that is strongly 
orbital dependent. Aex{dx2-y2) f« 1.3 eV is an unusually 
large value (LSDA exchange forces usually are less than 
1 eV//is), reflecting a weakly hybridized (and therefore 
more confined) dx2_y2 orbital. A^x is only about 0.6 eV 
in the center of the d bands and decreases to 0.4 eV at 
the top. The calculated gap in the majority bands is 0.7 
eV, and the gap between occupied majority and unoccu- 
pied minority bands is of the order of 0.1 eV. Correlations 
effects can only increase these gap values, perhaps sub- 
stantially. 





FIG. 3. Side view of the VO5 "pyramid." The V ion ac- 
tually lies above the center of mass of the five O ions. V-O 
bond lengths {A) are shown. The three V-01 lengths are 
inequivalent but nearly equal. 

We have applied the full potential linearized aug- 
mented plane wave (FLAPW) method ||l2|,^ used in 
previous studies of magnetic transition metal oxide com- 
pounds. The sphere radii used in fixing the FLAPW 
basis were chosen to be 2.0, 2.0, 1.60, 1.45, and 1.60 a.u. 
for Ca, V, 01, 02, and 03, respectively. Local orbitals 
were added to the basis set for extra flexibility and to al- 
low semicore states to be treated within the same energy 
window as the band states. For the V atoms, s, p, and 
d local orbitals were used, while for Ca and the O atoms 
s and p local orbitals were added. The plane wave cutoff 
corresponded to an energy of 19 Ry. The total basis set 
size is 3350 LAPWs and local orbitals. 

An initial calculation without spin polarization found 
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FIG. 4. Majority (left) and minority (right) bands along 
high symmetry directions for FM aligned CaV4 09. The ma- 
jority d^2_y2 is disconnected from other d states, leading to 
an insulating result. 

The spin density for this FM state, shown in Fig. 5, is 
clearly that due to occupation of a single spin-orbital of 
V dx2-y2 character, where x and y refer to the lines of V 
ions. Due to the low V ion site symmetry the lobes do 
not point precisely in the direction of neighboring V ions 
nor lie exactly in the crystallographic x — y plane, but 
these 'misalignments' are only a few degrees. These spin- 
orbitals are non-bonding with respect to the neighboring 
O ions, a point we return to below. This occupied orbital 
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is not of the type {dxy, dy^, dxz) anticipated in previous 
work, [^,|| or used by Marini and Khomskii for orbital 
ordering {d^z or dyz). Note that this orientation is 
determined by the crystal field, and is unrelated (at least 
for this narrow band system) to the FM order. 

The chosen AFM order was of the Neel type: each V 
spin is antiparallel to its two neighbors on a plaquette 
and to the neighbor on the next plaquette. This type 
of order breaks inversion symmetry, and all V ions on 
one side of the (idealized V-0) plane have the same spin 
direction. This ordered state is essentially degenerate 
with the FM alignment. The resulting state is also 
insulating, and as expected the occupied bandwidth is 
smaller (20%) than for FM alignment. 




FIG. 5. Top view of an isocontour of the spin density of 
FM aligned CaV409 (truncated at the cell boundaries), il- 
lustrating the occupied d^2_y2 V*"*" spin orbital. The white 
spheres indicate oxygen sites, small dark spheres indicate Ca 
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y2 transfer integrals ti 



are thought to be mediated via virtual hopping processes 
through the intervening O ions, with amplitudes tdpa- and 
tdp-K- In the idealized VO parent layer the O ions lie at 
45° angles to the V-V nn direction. In terms of hop- 
ping processes, the nnn V ion is no more distant that the 
nn V ion and therefore should have comparable hopping 
amplitude, hence exchange coupling J. According to the 
Goodenough-Kanamori (GK) rules, the hopping am- 
plitudes depend strongly both on the symmetry of the d 
and p orbitals involved and on the V-O-V angles. The 
first point to note is that this spin orbital is orthogonal 
(except for the « 5° rotation) to the O Pa orbital; hence 
the coupling must proceed through the orbital. This 
coupling involves only the tdpTr , which may be a few times 



smaller than tdpa- The Pt^ orbital couples to the d^2_y2 
orbital on the nnn V (again via tdp-,r)- The orbital 
however is orthogonal to the d 
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orbital on the nn V 
(with respect to which it is a Pa- orbital) . Thus the usual 
second order nn V superexchange coupling vanishes for 
the idealized layer. 

Due to the symmetry lowering distortions, this van- 
ishing is incomplete, and the surviving coupling may be 
larger than third order effects involving polarization of 
the O ion. pda hopping is allowed as tdpasincj) and 
pdvr hopping is reduced to tdp^cos4>. The couplings be- 
come 



Jnn OC 2 tdpjT C0S4> tdpa sin4>/ {Sd ~ Ep), 



(1) 

(2) 



where the factor of two arises from the two paths, and 
ed — Sp ~ 4-5 eV is the energy separation between V d 
and O p states. Since sint^ « 0.1, an rough estimation 

gives Jnnn/Jnn « \{t dpir / tdpa) COt (j) ~ \ ■ \ ■ ^ 2. (A 

rotation f2 necessary to rotate the dx2-y2 orbital into the 
plane of the V-0 bond affects both couplings proportion- 
ately, and has been neglected in this discussion.) 

This already rough estimate will be altered by the 
structural corrugations, nn V ions, which differ in z 
coordinate by 1.25 A, are connected along two V-O-V 
paths, whose contributions (presumably) add. For nnn 
V ions, which are at the same height, there is a single 
path. Within a plaquette, the angles for nn V ions are 
100° (V-Ol-V) and 95° (V-03-V), and is 144° (V-03- 
V) across the diagonal. Between plaquettes, nn V are 
connected by two identical 100° angles (V-Ol-V), while 
the diagonal nnn angle is 130° (V-Ol-V). These angular 
variations, coupled with the fact that the O ions are not 
coplanar with the lobes of the dj.2_j,2 spin orbital and 
that direct V-V exchange may not be negligible, make it 
very difficult to estimate realistic net exchange couplings. 

The direct V-V separations are 3.00±0.0lA, which can 
be compared to the bondlength in V metal of 2.73 A. 
This difference of less than 10% suggests that direct V-V 
hopping might be appreciable. The less extended na- 
ture of the V''^"'" orbitals compared to neutral V orbitals 
will reduce the overlap, however. Since nn V ions do 
not lie in the plane of the d^^-y^ orbital, the orbital 
must be rotated by 0=sin~^( j||)=25° on both atoms 
before the overlap can be expressed in terms of the usual 
tdda,tddir,tdds hopping amplitudes. Neglecting the lat- 
ter two, which at this distance should be much smaller 
than tdda, the effective hopping between the non-aligned 
dx2_y2 spin-orbitals will go as tdda cos^ 9, a reduction by 
cos^ e «0.8. 

It is not a straightforward matter to determine the 
relative importance of direct V-V hopping and V-O-V 
mixing that gives the dispersion pictured in Fig. 4. The 
effect on the d bandwidth of the V-O-V coupling was 
probed by moving the 02 ions from -=0.939 to -=0.75. 
This distortion brings the 02 ion nearer one V but much 
farther from the other (2.6 A), essentially eliminated V- 
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01-V coupling of neighboring plaquettes. The resulting 
occupied d bandwidth of W(i=0.4 eV, compared to 1.3 
eV for the real structure, can be ascribed to direct V-V 
coupling. This result suggests that coupling through the 

01 ions is responsible for roughly | of the bandwidth, 
with direct V-V interactions providing the remaining i; 
i.e. direct V-V coupling is not negligible. This complex- 
ity arises specifically from the fact that the spin-orbital 
is dj.2_y2 character, which greatly hinders coupling to O 
ions and maximizes V-V interaction. Using a coordina- 
tion number Zc=3, the V-V bandwidth translates to a 
hopping amphtude tdda =Wd/(2zc cos^ 6')=80 meV. 

The hybridization leading to magnetic coupling is re- 
flected in the spin density distribution on O ions neigh- 
boring V. The primary feature, in both FM and AFM 
alignments studied, is an antiparallel polarization of the 
apical 02 by perhaps as much as 0.2 (the spin density 
attributable to a given ion is not precisely defined). The 

02 ion, however, is not involved in exchange coupling. 
The 01 ions have moments about i as large, aligned 
parallel, while the 03 ion's moment is even smaller (zero 
by symmetry for AFM case). Thus the 01 site should be 
more important in the exchange coupling process than is 
the 03 site. 

These results - the occupied spin orbital of dx2-y2 
character, J„„„ > J„„, Jv-v t^O, degenerate FM and 
AFM alignments - suggest a picture of coupled two 
metaplaquette system. Each metaplaquette is comprised 
solely of V ions either above the plane, or below the plane; 
the upper metaplaquette is highlighted in Fig. 1. The 
two metaplaquette systems, each topologically equivalent 
to the original, are coupled by at least three, possibly 
competing, exchange couplings J„„, J4„, and Jy-y, but 
their values will not be easy to ascertain. The degener- 
acy of the FM and AFM alignments, in which the two 
ferromagnetic metaplaquette systems are aligned and an- 
tialigned respectively, suggests relatively weak net cou- 
pling between the two metaplaquette systems. 

There is one particularly interesting difference between 
the original plaquette system and (either one of) the 
metaplaquette systems. As was noted in the introduc- 
tion, V-O-V couplings within a plaquette, or between 
neighboring plaquettes, are very similar from the struc- 
tural chemistry point of view. The metaplaquette sys- 
tem is different: intra- metaplaquette coupling proceeds 
through the 01 site, while inter-metaplaquette coupling 
proceeds through the 03 site. Not only is the V-03 dis- 
tance 4% longer than the V-01 distance (Fig. 3)result- 
ing in weaker coupling, but the environments are dif- 
ferent and the LSD calculation indicates substantially 
less hybridization through the 03 site, as reflected in 
the much lower induced 03 spin density (somewhat vis- 
ible in Fig. 3). This 01-03 site distinction provides 
a mechanism whereby the tendency to form singlets on 
each metaplaquette is not disrupted by equally strong 
inter-metaplaquette coupling, which was the case for the 
original plaquette model. 



To summarize, it has been found that consideration of 
the full low symmetry structure of CaV409 indicates that 
five unequal exchange couplings (Jnn, Jnnn, Jnn^ JLm 
from superexchange, and direct Jy-v) are present. The 
orbital character of the spin, together with the topology 
of the V-0 layer, suggest that a coupled two metaplaque- 
tte system may provide a realistic way to interpret the 
spin gap behavior. 
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